Effects of Axion-Photon Mixing on Gamma-Ray Spectra from Magnetized 

Astrophysical Sources 



Kathrin A. Hochmuth 1 and Giinter Sigl 2 

1 Max- Planck- Institut fur Physik (Werner-Heisenberg-Institut) , Fdhringer Ring 6, 80805 Miinchen, Germany 
" 2 APC * (AstroParticules et Cosmologie), 10, rue Alice Domon et Leonie Duquet, 75205 Paris Cedex 13, France 
and Institut dAstrophysique de Paris, UMR 7095 CNRS - Universite 
Pierre & Mane Curie, 98 bis boulevard Arago, F-75014 Paris, France 
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conditions under which 7— ray spectra from such sources would be significantly affected by axion- 
photon mixing. We show that, depending on strength and coherence of the magnetic field, axion 
couplings down to ~ (10 13 GeV) -1 can give rise to significant axion-photon conversions in the envi- 
ronment of accreting massive black holes. Resonances can occur between the axion mass term and 
the plasma frequency term as well as between the plasma frequency term and the vacuum Cotton- 
Mouton shift. Both resonances and non-resonant transitions could induce detectable features or 
even strong suppressions in finite energy intervals of 7— ray spectra from active galactic nuclei. Such 
effects can occur at keV to TeV energies for couplings that are currently allowed by all experimental 
constraints. 
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I. INTRODUCTION 

Axions and axion-like particles couple to two photons 
and can thus convert into on-shell real photons and vice- 
versa in the presence of magnetic fields. Astrophysical 
sources produce high energy 7— rays as secondary parti- 
cles of primary charged cosmic rays. Since these primary 
cosmic rays have to be accelerated in electromagnetic 
fields, the spectrum of any astrophysical 7— ray source 
can potentially be modified by axion-photon mixing. In 
the present paper we discuss in detail in which source 
environments and for which axion masses and coupling 
constants one can expect significant modifications of the 
observable 7— ray spectrum. As specific examples we con- 
sider the jets and central engines of active galactic nuclei 
(AGN) both of which are driven by magnetized accretion 
disks. 

We consider the Lagrangian for the coupling between 
a pseudo-scalar field a, called axion in the following, and 
the electromagnetic field strength 

C ia = --g ia F l j, v F' iV a = g ia E-'Ba, (1) 

where F^ v = ^Sf ll/pa F pa is the electromagnetic dual, E 
and B are the electric and magnetic field strengths, re- 
spectively, and g ia is the photon-axion coupling. Eq. ([1]) 
implies that a photon can convert into an axion in a mag- 
netic field B. In the present work we will consider the 
effect of such conversions on 7— ray fluxes emitted from 
magnetized sources for general pseudo-scalar masses m a 
and photon-axion couplings g ia - 
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Axions with couplings g ia < 10 _10 GcV~ 1 are con- 
sistent with all existing constraints for almost all axion 
masses [l| and we will use g ia = 1CP 11 GcV -1 as our 
benchmark value. As possible 7-ray sources we consider 
already detected sources, like Makarian 501 and Makar- 
ian 421. Sources of this type can have a model-dependent 
magnetic field strength between 0.01 G and several G. 
Adopting these values we find that a significant amount 
of 7-rays can be converted into axions, thereby leaving 
characteristic signatures like steps or gaps in the spec- 
trum. This conversion can take place due to normal os- 
cillations and due to resonance effects. We find that for 
typical parameters the effects can be sizable and are well 
within the detection region of the next generation of ex- 
periments. 

Our paper is organized as follows. In Sec. [Til we discuss 
in general how the 7-ray spectra can be modifed due to 
resonant and non-resonant oscillation effects. We discuss 
the prospects of detecting these effects in Sec. IIIII where 
we discuss specifically the case of central engines and jets 
of AGN. In Sec. IIVI we present our conclusions. We use 
natural units, h = c = k = 1, throughout the paper. 

II. CONDITIONS FOR MODIFICATION OF 
GAMMA-RAY SPECTRA 

The axion-photon conversion probability in a trans- 
verse magnetic field B t can be derived by using the lin- 
earized equation of motion following from Eq. ((TJ) for rel- 
ativistic axions. This equation can be written as 0] 

(A ± \ 

(E-id z -M) I A,, 1 =0, (2) 

where z is the direction of propagation, E is the photon 
energy and a the axion field. Aj_ and Au arc orthogonal 
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components of the photon field, where i = _L or || refer 
to the B t direction. The mixing matrix M. is 
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Here we have defined lu 2 ^ = Aira n e /m e as the plasma fre- 
quency for an electron density n e , where m e is the elec- 
tron mass and a = e 2 /(47r) the fine-structure constant 
with e the electron charge. Ar is the Faraday rotation 
term, which is dependent on the energy of the longitudi- 
nal component of the magnetic field and couples to the 
modes Au and A_\_. However, we want to consider only 
non-polarized sources, which renders this term negligi- 
ble. The vacuum Cotton-Mouton effect is represented 
by the Acm terms, which describe the birefringence of 
fluids in presence of a longitudinal magnetic field, with 
|Aq M — Aq M | oc B 2 . Note that the photon dispersion 
relation for polarization i is Ei(k) = k + Aj, correspond- 
ing to the refractive indices rii — 1 ~ Ai/E. For ax- 
ions one has, of course, E(k) ~ k + m 2 /(2k) and thus 
A a ~ml/(2E). 

Neglecting different photon polarization states and de- 
noting the resulting photon state by A, we are left with 
a simple two-component mixing problem, 



E - id z — 



Api(z) + Acm(z) Afl(z) 
A B (z) A a 



= 0, 
(5) 

where we have indicated the terms that are location de- 
pendent. After diagonalization of the mixing matrix we 
obtain the solution 
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where E is the photon energy and B cr = m 2 /e ~ 4.41 x 
10 13 G is the critical magnetic field strength. Since we 
are only interested in orders of magnitude, we neglect dif- 
ferent photon polarization states whose vacuum Cotton- 
Mouton term have slightly different prc-factors, and we 
have averaged over directions, (B 2 ) ~ B 2 /3. Further, in 
Eq. © we have used the abbreviations B m Q = (-Bt/mG), 
.9ii = 9ya x 10 11 GeV, and ra Me v = (m a /fieV). 

Only terms to lowest non-trivial order in B are taken 
into account in Eq. ©. The term to next-higher or- 
der in B in Acm is suppressed by a factor ~ x 2 = 
[(-E/TOeX-B/Bcr)] 2 @ ■ This is negligible for energies 
£7 < 2 x 10 16 (kG/B) eV. Furthermore, the magnetic field 
influences the phase space distribution of the electrons 
which modifies the plasma term to ~ (E/2)ui 2 l /(E 2 — 

w c) @, where u> c — eB/m e = (B j B cr )m e is the elec- 
tron cyclotron frequency in a non-relativistic plasma. 
The latter is thus negligible for energies E > lo c ~ 
10~ 5 (73/kG) eV. Furthermore, inelastic processes such 
as pair production and photon splitting in the magnetic 
field are suppressed as exp(— l/x) an d thus negligible for 
X C 1 @|. Since we are mostly concerned with magnetic 
fields B < 10 4 G and photon energies E < PcV, both in- 
elastic processes and higher order corrections to Eq. © 
can be neglected. 

For propagation over TV coherence domains the total 
conversion probability of photons into axions can be de- 
rived as 0] 



[1 - cxp (-3iVP /2)] , 



(10) 



Within a domain of linear size s and roughly constant 
plasma density and magnetic field, the probability of a 
transition from a photon to an axion is 0, [H, 0| 



which saturates to P 7 ^ a ~ | for NP 3> 1. 
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where the oscillation wave number is 
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In Eq. ([5]), the plasma contribution A p i, the vacuum 
Cotton-Mouton term Acm, the axion mass term A a , and 
the off-diagonal mixing term Ag are given by 



A a = — 
2E 



2.5xlO- 20 m 2 eV 



TeV\ _ x 
— — cm 



A. Resonances 

In the following, let us denote the coherence length of 
the plasma and magnetic field by A with A = A/pc. The 
coherence length is of course always smaller than the size 
of the system. 

MSW type resonances can occur in Eq. ([9]) when A p i 
becomes comparable to A Q or the modulus of Acm and 
the other contribution is negligible. In order for the res- 
onance to lead to efficient conversion, it has to be adia- 

< A 2 B at the resonance, 
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batic, requiring 2n 
where the prime denotes the derivative with respect to 
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the distance along the propagation direction. In the fol- 
lowing we will estimate this derivative by the coherence 
length, i.e. ' — > 1/A. 

One has A a = A p i at the electron density 



n e ~ 7.1 x 10 8 m 2 cV cm 



(11) 



This leads to an adiabatic resonance at energies satisfying 



E > 20g- 1 2 ml cV \- c 1 B^ i T C V. 
E < 4.4m^v-B-^TeV. 



(12) 



In Eq. (|T2j) , the first condition results from the adiabatic- 
ity requirement with A' osc ~ A p i/A and the second from 
Acm ^ A a . As a result, strong 7— ray flux suppression 
is expected at energies satisfying Eq. (TT2")) . 
One has Acm = — A p i at the energy 



E ~ 5.2 



VlO 3 cm- 3 / 



1/2 



(13) 



This leads to an adiabatic resonance provided that at this 
energy A a < |A p i| and that the adiabaticity condition is 
fulfilled with A' osc ~ Acm/ A. This translates into 



(14) 



n e > 7.1 x 10 8 m 2 cV cm 3 
E < 0.95 A pc TcV, 



respectively. Note that the first condition in Eq. (TT4"|) 
together with Eq. flT3J) implies E > 4.4 m^ cV B~^ G TeV, 
such that this type of resonance starts to be relevant 
at energies just above where the first type of resonance 
ceases. This is understood since for the first resonance 
we had |Acm| iS A a , whereas for the second A a < A p i = 
|A CM |. 

Eq. (fi"2"|) and the second inequality in Eq. (| 14|) com- 
bined thus lead to resonances in the energy range 

E > 20 5 riXovA pc 1 S mG TeV, 

E < 0.95.9?! A pc TeV, (15) 

provided the following conditions are satisfied: 

1/2 
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Here the first condition comes from the requirement that 
Eq. (|15p represents a finite energy range, and the second 
condition results from Eq. (jlip and the first inequality in 
Eq. (|14[) . For a given coupling gn resonances occur only 
for coherence lengths (and thus system sizes) satisfying 



d > A > 1.4 X 10 19 g^m^vB^ cm 



(17) 



Inserting Eq. (fT7|) into Eq. (|15p thus also leads to the 
inequalities 



E 1 ■ < 4 4 



"fyeV £> mG 



TeV < EL, 



(18) 



for the minimal and maximal photon energy E r min and 
E* x , respectively, at which resonances occur. 



B. Non-resonant Oscillations 

Let us denote the coherence length of the magnetic 
field by A and the propagation length with d, so that 
we have N ~ d/X > 1 domains. For A osc A < 1, we 
have from Eq. (JT]) Po ~ (A^A) 2 . In the opposite limit, 
A osc A > 1, we have Po ~ 2(Ab/A osc ) 2 . Altogether, with 
Eq. (jTUJ) this yields for the photon survival probability 



p 



1 - P 



7— >aj 



exp 



min [AldA^A^As/Aosc) 2 



(19) 

For given axion parameters, the spectra predicted within 
scenarios of 7— ray sources will be modified by this factor. 
Eq. (TlT)|) can only be of order unity and thus lead to ob- 
servable effects, if both A 2 B dX > 1 and 2(A B /A osc ) 2 A > 
1. With Eq. the first condition yields 



(dX) 1/2 > 5.7 x 10 19 g^ 1 B~ G cm . (20) 
This results in the two conditions 



and 



d>5.7xl0 19 g^ 1 B- G cm 1 
■d^ 1/2 
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1 < N 1 ' 2 = ( - ) < 0.054 d pc g n B mG , (22) 



or equivalcntly 



A 



-> 320 3 r i V B mGPC 



(23) 



where d pc = (rf/pc). 

The second condition from Eq. (TT5]) leads to the 
three separate inequalities As <; A a A -1 / 2 , As ^ 
I A p i I iV" 1 / 2 , and A B > A CM iV~ 1/2 . Substituting Eq. © 
results in 

E > 15 g^m 2 JV-^B-^TeV, 



E > 21 



10 3 < 



g^ N- 1 / 2 B- G MeV , (24) 



E < 1.3ff U JV 1 / 2 B-^'lteV. 

Note that the scaling of these energies with A are different 
from the scaling of the resonance energies Eq. (fTS"]) : Large 
coherence lengths, or small N, favor a broad resonance 
energy range, whereas small coherence lengths will tend 
to lead to non-resonant transitions, as long as Eq. ([22]) is 
satisfied. The absolute minimal and maximal photon en- 
ergies at which non-resonant oscillations occur are given 
by setting N = 1 in Eq. l[2"4"]). 

We can now eliminate N or equivalently A from 
Eq. (HD by using Eq. 



E > 278.gr 1 2 m 2 cV d- 1 p- 2 3 TeV, 



E > 389 
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(25) 



E < 0.07 5^ d pc TeV . 
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For Eqs. (|24|) and (f25|) to be satisfied for a finite energy 
range the following conditions are implied: 



9u 



i > 3.4m„ oV N" 1 / 2 > 0.13 



m^ev 



(26) 

,ttpci>mG / 

< 6.2 x 10 7 .gi! N cm~ 3 < 1.8 x 10 5 g^c^mC cm -3 , 



1/2 



where in the second expressions we have substituted 
Eq. (|2"2"|) . The condition on <?n is very similar to the 
resonance condition Eq. (fTo| . 

Coupling and mass of the QCD axion roughly satisfy 
<7ii ~ 10 m^ e V) and thus experimental limits imply 
m a ^ 1 cV. The first condition in Eq. (|26p then implies 



N > 1.2 X 10 11 , 



(27) 



rf P cS mG > 18 5^ - 1.8 x 10 6 m-V 



where in the second line we have used Eq. (f2"Tj) . In 2005, 
the PVLAS Q experiment has seen indications for axion- 
photon mixing with 



m a ~ 1.3 mcV 

ff7a ~ 3 x 10~ 6 GeV -1 , 



(28) 



which since recently, however, is considered to proba- 
bly have been an experimental effect @. In any case, 
the "PVLAS axion" would fulfill the first condition in 
Eq. (EH). 



C. Other Conditions 

A further condition comes from the requirement that 
the length scale over which the emission is created can- 
not be larger than the variability time scale of the source. 
The size of the emission region is determined by the 
length scale over which the plasma becomes transpar- 
ent to 77 pair production. Emissions of active and ra- 
dio galaxies often vary on times scales of days or even 
hours [{J [TTJ], corresponding to scales < 10 14 ~ 15 cm. If 
these are the same length scales over which significant 
photon-axion conversion occurs, then one also has the 
condition A < d < 10 14-15 cm, except in case of relativis- 
tic beaming. 



III. PROSPECTS FOR SPECIFIC SOURCES 

A lower limit on the source magnetic field times the 
propagation length can actually be obtained by requiring 
that it accelerates cosmic rays of charge Ze up to the 
maximal energy E CI observed from that source: 

d > 1^3.3* 10^(1^. on, (29) 

This implies that sources that accelerate particles up to 
E CI >1.7x 10 19 Zg^ 1 cV (30) 



fulfill the condition Eq. (|2Tj) independently of the mag- 
netic field strength. The coherence length of the fields 
depends on the acceleration model and has to satisfy the 
additional condition Eq. (j2li|) . Sources of ultra- high en- 
ergy cosmic rays with energies E'er > 10 18 eV are thus 
promising objects that may exhibit axion-photon mixing 
induced modifications of their photon spectra for cou- 
plingsgn > 1.7(Z10 19 eV/£ cr ). 

Whereas this condition for the occurrence of photon- 
axion oscillations does not depend on size or magnetic 
field strength of the accelerator, the energies at which 
such oscillations could modify photon spectra do depend 
on these parameters: Eqs. (|15[) and (|25[) show that the 
minimal and maximal energies at which significant oscil- 
lation effects can occur, are given by 



E n 



,'Z10 18 eV\ 
IQSnXcV^mG f ) TeV 



183ii 2 m^ eV d pc 



^cr 

Z10 18 eV 

E rr 



TeV 



0.95 ^dpc TeV ; 



(31) 



where we have used the relation Eq. (j29l) . 

Among the extragalactic objects which have been seen 
in very high energy 7— rays and where photon-axion con- 
version could occur are Markarian 42 1 [ill . Markarian 
501 [HI, [T3L the blazar 1ES 1101-232 (Tjand the vari- 
able core Q and flaring knots Y\M of M87. Markarian 
421 is also variable in 7— rays[10j- Recent reviews on 
leptonic and hadronic models of blazar emission can be 
found in Refs. [HI, [l6j]. The spectrum of a typical AGN 
has a double-peaked power spectrum. In leptonic mod- 
els, the high energy peak is caused by inverse Compton 
scattering of accelerated electrons on the ambient photon 
field, whereas the low energy peak is due to synchrotron 
emission of these same electrons. The energies and rela- 
tive power flux in these peaks thus contains information 
on the magnetic fields in the emission region. The esti- 
mated magnetic field strengths range between ~ 10 mG 
and ^ 10 G over scales, depending on the models used, 
typically of order d > 10~ 3 pc. This would also suggest 
possible cosmic ray acceleration up to ~ 10 18 eV. For ax- 
ion parameters gn ~ 1, m^y ~ 1, Eq. (f3"Tj) then suggests 
photon-axion oscillation effects_at GeV energies, consis- 
tent with the findings of Ref. 



17] 



For the galactic center fields of order 10 G and up to 
10 4 G have been discussed in Ref. [3, UH. The galactic 
center is sometimes thought to accelerate cosmic rays up 
to - 10 18 eV [2(|. Eq. J31]) would then imply significant 
non-resonant axion-photon oscillations down to MeV en- 
ergies. 



A. Former Work 

The original indications from PVLAS for axion-likc 
particles with rather strong couplings, Eq. (f2"8"l) . although 
in contradiction with present astrophysical constraints, 
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has motivated the study of several possible astrophysi- 
cal effects. Ref. HH, H3, considered variable 7— ray light 
curves from double pulsars as signatures for the PVLAS 
axion-photon mixing parameters. In this scenario, a sur- 
face pulsar magnetic field B ~ 10 12 G yields B ~ 10 4 G, 
TV ~ l in the accretion region between the pulsars 
where the 7— rays are produced, and the Goldrcich- Julian 
plasma density is n e ~ 10 3 cm~ 3 . Eq. (j2"4")l implies a 
modulation for MeV < _E < 100 GcV, as obtained in 
Ref. [21], [22j ■ Note that for less coherent fields between 
the two compact stars, TV > 1, the modulation would 
extend down to lower energies. 

Ref. Q considered the modification of the galactic cen- 
ter 7— ray flux by PVLAS axions, for B ~ 71G, A ~ 
0.01 pc, d ~ lOkpc, thus TV - 10 6 . Thus Eq. results 
in a modification at energies lOTeV < E < 10 12 TeV. 
Note that the coherence length of the Galactic mag- 
netic field is not very well known [22j. For A > pc, 
the effect would disappear for all observable energies 
E < 100 TeV. Measurement of Faraday rotation from 
pulsar pairs suggest that the cell size of the random com- 
ponent of the gala ctic magnetic field may indeed be of 
order 50 pc [Sill- 

No resonances are expected in these two scenarios due 
to the second condition in Eq. (fl6|) . unless m < 10 -8 eV. 

Ref. [2(| constrained axion-photon mixing by using 
quasar spectra mostly in the optical, considering mixing 
in intergalactic fields of strength B < 10~ 9 G only. In 
this case such spectra are only sensitive to axion masses 
m a < 10~ 15 cV, consistent with Eqs. (JT5J) and p5]). For 
larger axion masses m ~ 10~ 10 cV, such conversions in 
intergalactic magnetic fields would show up at TeV ener- 
gies, as discussed in Ref. [27 1. 

Ref. [13] discussed the modification of 7— ray spec- 
tra from AGNs by axion-photon mixing, using magnetic 
fields and length scales according to the Hillas criterion 
for accelerating particles to the pertinent energies. How- 
ever, they did not consider resonances. 

Note that emission from the surface of a neutron 
star or magnctar with d ~ 10 6 cm, B ~ 10 12 -10 15 G, 
n e ~ 10 34 cm -3 according to the second condition in 
Eq. (|26p could not give significant non-resonant transi- 
tions for <7n < 10 5 . Furthermore, according to Eq. (fTS"]) . 
resonances can only occur at very low energies E < 
0.3 ^ eV. 



B. The central engine of AGNs 

Black hole accretion launches jets by magnetohydro- 
dynamical processes. For spherical Bondi-accretion, the 
accretion rate is M ~ 4irr n e m p P where m p is the proton 
mass and (3 ~ Poifs/f) 1 ^ 2 ar e the density and velocity 
of the accretion flow, respectively. Here, (3q ~ 0.1 and 
r s = 2G N Af ~ 2.96 x 10 14 Af 9 cm is the Schwarzschild 
radius for a black hole of mass M = Mg M Q . In- 
troducing the Eddington luminosity, LEdd(M) — 1.3 x 
10 47 Mg ergs -1 , the accretion rate can be written as 



M = fEddLEdd(M)/rj, where 77 is the efficiency with 
which accretion is converted into electromagnetic radi- 
ation, and /Edd = L\y \/ L~Edd{M) is the ratio of the bolo- 
metric luminosity to the Eddington luminosity. AGNs 
have duty cycles of a few percent and in their active pe- 
riods have 77 ~ / Edd ~ 0.1. 

From this we can estimate the plasma density in the 
accretion flow by n e ~ M/(47rr 2 m p /3) which close to the 
Schwarzschild radius yields 



2.6 x 10 



9 ( /] 



Edd 



M 



rs 



-3/2 



(32) 



If the magnetic fields arc roughly in equipartition 
with the accretion flow at radius r, one has B 2 /(8w) ~ 
n e m p p 2 . This results in [28| 
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where in the second expression the accretion rate has 
been written in terms of the Eddington luminosity. Note 
that these are very rough order of magnitude estimates 
which, however, will be sufficient for our purpose. 

If accretion is non-spherical, but instead occurs in a 
disk, as is in general the case for very luminous AGNs, 
for a given accretion rate, the plasma density and mag- 
netic field strengths tend to be larger than the estimates 
Eqs. (|3"2"|) and (j3"3"|) . In these and all subsequent formu- 
las, this is essentially mimicked by an electromagnetic 
efficiency r\ smaller by a factor of roughly the height of 
the disk relative to its radius. This will be very roughly 
the Lorentz factor of the resulting jet. The geometry is, 
however, more complicated in this case and photon-axion 
oscillations may be suppressed perpendicular to the disk. 
The strongest effects may then only occur in directions 
of the plane of the disk. As we will sec below, both the 
minimal Eddington ratio / E dd and the energy for which 
oscillation effects occur are proportional to 77, whereas 
other parameters are largely insensitive to r\. 

In any case, more detailed theoretical models of mag- 
netized accretion disks lead to magnetic fields of or- 
der 10 4 G in the disk and - 10 2 G above the disk 
Such values are consistent with the rough estimates in 
Eq. (f33|) . In addition, supermassivc black holes may have 
a dipole magnetic field with a maximal field strength of 
- 2 x 10 10 Aig" 1 G at the Schwarzschild radius [U. 

But how coherent are these fields? We note that they 
may be made coherent by dynamo effects and/or mag- 
netorotational instability (MRI) Jet launching by MHD 
effects may indeed require coherent, very non-thermal 
fields. In this case the fields would be coherent on roughly 
the Schwarzschild radius. 

Magnetic fields of order 10 G in the vicinity of the 
radio galaxy M87 have been discussed in Ref. [31j|. This 
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value is consistent with the estimate Eq. (|33j) if M87 is 
in a low accretion state, / E dd ~ 10~ 6 , which is the case 
according to present measurements and theories [3ll |. 

The variability constraint is fulfilled automatically 
since AGNs cannot be variable on time scales smaller 
than the Schwarzschild radius. 



1. Resonances 

Comparing the second condition in Eq. (fl6 
Eq. (|3"!?|) implies the condition 

M<3.7xl0 9 < v M 



with 



(34) 



for a resonance to occur above the Schwarzschild radius. 

Using Eq. (|3~3"1) and assuming A « r, Eq. (fTo*)) shows 
that the widest energy range is achieved for the smallest 
radii. Further, the first condition in Eq. (|16|) implies the 
least stringent condition on gn for the smallest radii. We 
thus get 

SU > 8.8 x 10- 2 m% M 9 ^ (-5-) ^ . (35) 



Thus, when the above conditions Eq. (|34|) and Eq. (|35|) 
arc fulfilled, a strong 7— ray flux suppression occurs for 
an energy interval given by substituting the field strength 
Eq. (|3"3"1) close to the Schwarzschild radius into Eq. (|15|) . 



E > 6xl0 3 5n 2 m 2 oV iV (— 
£ < 90 g^N- 1 M 9 MeV, 



Edd 



eV. 



(36) 



where iV = A/rg is the number of coherent domains 
within a Schwarzschild range. Note that the lower limit 
does not depend on the black hole mass. In Fig. [1] we 
show the energy range Eq. (|3"6"|) for fields coherent over 
the Schwarzschild radius, N = 1, for an AGN of the type 
of Mrk421 or Mrk501, with M ~ 10 6 M Q , / Ed d ~ 10~ 3 . 
For couplings approaching the experimental upper limit, 
9a-y iS 10~ 10 GeV" 1 , resonances can occur down to ener- 
gies ~ 600 (m/0.1/ieV) 2 (/Bdd/10" 3 ) _1 eV which can ex- 
tend into the keV regime. 

The condition on the coherence scale Eq. (fT7j) reads 



/Edd > 6.1 xlO- 5 (^-)g u 4 ml cV M g 



N- 1 = —> 7.8 x 10- 3 .g 1 - 1 2 m /ie vM 9 
rs 



-1/2 



Pofi 



Edd 



(37) 

1/2 



where the first condition follows from the second observ- 
ing that N > 1. 



2. Non-resonant Oscillations 



Substituting d = r 
Eq. (HI]) and in Eq. 



rg in the conditions Eq. (J^J) and 
results in requirements similar 



> 




FIG. 1: Maximal and minimal energies Eq. ()36)l for which 
axion-photon resonances occur, for M ~ 1O 6 M , / E dd ~ 
10~ 3 , /So = n = 0.1, A = rs. Whenever E m in < E msx , 
adiabatic resonances occur. Note that in this case also the 
conditions Eqs. ([34]), ([35j, and ([371) are fulfilled. 



to Eq. §7 



/ Edd > 1.1 X 10~ 3 (j-jg^Mg 1 , 

N- 1 = A>l.lxl -3 5ri 2 M -l 

PO/Edd 



(38) 



'/ 



We remark that the first constraint is easily fulfilled for 
luminous AGNs which have /Edd ^ 10 -2 , especially for 
relatively large axion-photon coupling. The central black 
hole of our Galaxy has / E dd ~ 10~ 8 . 

Note that the condition on the coherence length in 
Eq. (|3"8)) tends to be less stringent than Eq. ([3"T]) . con- 
sistent with the fact that more coherent fields tend to 
lead to resonances. 

Using the plasma density Eq. (|3"12|) and again substi- 
tuting d = r = rs and Eq. (|33p . the two conditions in 
Eq. (J2H) translate into 



gu > 5.5xl0- 3 my c 2 v Af 9 1/4 
gn > 0A6f3~ 1/2 M g 1/2 . 



A) /Edd 



1/4 



(39) 



This implies that emission from massive black holes can 
be sensitive to rather small coupling constants. The 
range of coupling constants given by Eq. (|39p is shown in 

Fig. m 

Thus, when the above conditions Eqs. (|38|) and 
Eqs. (|39[) arc fulfilled, a 7— ray flux suppression by about 
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1000 



100 r 




FIG. 2: Range of coupling constants and axion masses 
Eq. (|39p for which non-resonant axion-photon oscillations can 
influence the 7— ray spectra from AGNs with /Edd = 10~ 3 . 



a factor ~ 2/3 occurs for energies 



1/2 



M 9 1/2 MeV 

PO/Edd 



> 0.083 g^i 



2 _2 



Edd 



MeV, 



(40) 



E < 0.22 g^N 

< 6.6 M 9 MeV, 



1 f _J? 



Edd 



1/2 



M Q 1/2 MeV 



where for the second expressions we have used the con- 
straint on the number of domains in the second line of 
Eq. (|38| . The third condition that can be obtained from 
Eq. |24|) is less stringent than the one given. 

Note that the energies Eq. (|40p tend to be in the hard 
X-ray to soft 7— ray range except for very strong cou- 
plings of the PVLAS type for which these energies can 
extend down to the optical range. Such photons are not 
absorbed by pair production and there is plenty of data 
at such energies that could be searched for the spectral 
features discussed here. 



C. AGN jets and hot-spots 

AGN jets and hot-spots seem to be likely candidates 
to fulfill the requirements for observable photon-axion 
conversion, as the condition Eq. (|2Tj) is satisfied for all 
couplings gn > 10~ 2 in such environments since even 
the transverse dimension of the jets are of order kpc and 
fields are at least of order 100 ^G [H, SI, [H| . The first 



condition in Eq. (|26|) shows that this allows significant 

1 /2 

effects down to couplings g\\ > 0.013 m^y, provided the 
coherence scale A ~ pc. Furthermore, the plasma density 
in jets is assumed as n e < 10 5 cm~ 3 , about a factor 100 
higher than the ambient density of an average galaxy, 
as expected for termination shocks. Thus, the second 
condition in Eq. (f2"6")) is also satisfied, unless gn < 0.1. 
According to Eq. ([23)1 , the resulting features would show 
up at TeV energies. However, due to Eq. 1]14[) . resonances 
would only occur for m a < 1.2 x 10~ 8 eV. 

In the proton synchrotron model, field strengths of up 
to 10 mG have been discussed in Ref. [35j]. For axions 
with gu ~ 1, m a ~ /LieV Eqs. (|24|) and (|25j) thus imply 
possible effects down to ~ GeV energies, if A ~ 0.1 pc. 



1000 r-? 



a 



T3 




0.001 0.01 



10 100 
E [GeV] 



1000 10 4 10 5 



FIG. 3: Modification of the 7— ray spectrum of an AGN of 
mass M = 10 9 M Q , for B = 0.5 G over 2 x 10 6 Schwarzschild 
radii, or ~ 200 pc with coherence length A equal to 10 times 
the Schwarzschild radius, or ~ 10 -3 pc. The assumed in- 
^-2.25 ax j on p aram eters are 



jection spectrum is oc 



10" 



GeV" 



P 7 _ 7 is given by Eq. 



m a = 1/ieV. The modification factor 
(|19|) . Spectral modification by pair pro- 



duction has not been taken into account and is negligible for 
redshifts z < 0.03 

As an example for the effect of non-resonant oscilla- 
tions we show in Fig. [3] the modification of the 7— ray 
spectrum of a quasar around the maximal energy given 
by Eq. (f2"4"]) . The parameters used are typical for AGNs 
with central black hole masses around 10 9 M Q . The mag- 
netic field is assumed of order a Gauss over length scales 
10 6 times the Schwarzschild radius. Such fields can oc- 
cur in models of kpc scale jets which emit MeV 7— rays 
produced by inverse Compton scattering of accelerated 
electrons on low energy synchrotron photons and exter- 
nal photons [lllHIl; an d in ultra-compact jets [13] • How- 
ever, variability over scales of months or less would imply 
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that the length scales over which significant photon-axion 
conversion occurs would be much larger than the size of 
the region where most of the emission is produced. This 
region would have to be close to the black hole, except 
if the jet is oriented toward the observer such that rela- 
tivistic beaming compresses the time scales. 

IV. CONCLUSIONS 

We have investigated the possibilities for successful 
photon-axion conversion in 7-ray sources such as the discs 
and jets of AGNs. We have derived conditions for reso- 
nant and non-resonant oscillations and find that a signifi- 
cant conversion is possible for standard axion parameters 
with coupling gn ~ 1 and axion mass m^y ~ 1, which 
are allowed by present experimental and astrophysical 
constraints. Of course the efficiency of conversion de- 
pends on the strength of the magnetic field in the vicinity 
of the black hole. Values between ~ 0.01 G and several 
G are under discussion. Resonant effects, leading to gaps 
in the observed spectrum, are observable for plasma den- 
sities n e > 7.1 x 10 s m 2 ^ cV cm -3 and energies between ~ 
keV and MeV for emissions from the central engines of 
active galactic nuclei. Such resonances can occur between 
the axion mass term and the plasma frequency term as 
well as between the plasma frequency term and the vac- 



uum Cotton-Mouton shift. Assuming a limited region of 
production and conversion of order the gravitational ra- 
dius we find that non-resonant effects can occur again in 
the keV and MeV energy range, provided that also the 
magnetic field is coherent over this distance. AGN jets 
and hot-spots should provide an even more interesting 
site with field strength of order /xG-mG and jet sizes of 
order kpc. In these scenarios we find that a significant 
conversion is possible for MeV-TeV energies. Our find- 
ings should be testable in high energy 7-ray experiments. 
Additionally, though all our considerations encompassed 
normal axion parameters, our limits can be applied to 
axion-like particles like the PVLAS axion. In this case ef- 
fects should even be observable down to the optical scale 
and thus be already testable. 
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